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[i] In this study, Moderate Resolution Imaging Spectroradiometer (MODIS) fine mode 
fraction and Multi-angle Imaging SpectroRadiometer (MISR) nonspherical fraction data are 
used to derive dust and smoke aerosol optical thickness (idust an d x smoke ) over the tropical 
Atlantic in a complementary way: due to its wider swath, MODIS has 3M times greater 
sampling than MISR, but MISR dust discrimination is based on particle shape retrievals, 
whereas an empirical scheme is used for MODIS. MODIS and MISR show very similar dust 
and smoke winter climatologies. x dust is the dominant aerosol component over the tropical 
Atlantic, accounting for 40-70% of the total aerosol optical thickness (AOT), whereas x smoke 
is significantly smaller than x dust . The consistency and high correlation between these 
climatologies and their daily variations lends confidence to their use for investigating the 
relative dust and smoke contributions to the total AOT variation associated with the 
Madden-Julian Oscillation (MJO). The temporal evolution and spatial patterns of the x dus 
anomalies associated with the MJO are consistent between MODIS and MISR: the magnitude 
of MJO-realted x dust anomalies is comparable to or even larger than that of the total x, while 
the x smoke anomaly represents about 15% compared to the total, which is quite different from 
their relative magnitudes to the total x on the climatological time scale. This suggests that dust 
and smoke are not influenced by the MJO in the same way. Based on correlation analysis, dust 
is strongly influenced by the MJO-modulated trade wind and precipitation anomalies, and can 
last as long as one MJO phase, whereas smoke is less affected. 
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1. Introduction 

[ 2 ] The Madden-Julian Oscillation (MJO) [ Madden and 
Julian, 1971, 1972] is the dominant form of intraseasonal 
(30-90 day) variability in the tropical atmosphere. It is charac- 
terized by slow (~5ms _1 ) eastward-propagating, large-scale 
oscillations in the tropical deep convection over the equatorial 
Indian Ocean and western Pacific during boreal winter 
(November- April) [Lau and Waliser, 2005; Zhang, 2005], 
Recently, there is strong emerging interest in the impacts of 
the MJO on atmospheric composition [Tian and Waliser, 

2011] , such as aerosol [Tian et al., 2008, 2011; Reid et al., 

2012] , ozone [Tian et al., 2007; Weare, 2010; Li et al.. 
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201 1], carbon dioxide [Li et al., 2010], and carbon monoxide 
[Wong and Dessler, 2007]. 

[ 3 ] Tian et al. [2008] first examined the aerosol variability 
related to the MJO using global aerosol products from mul- 
tiple sensors on various satellite platforms. That study 
revealed large intraseasonal variations in the satellite-derived 
aerosol products over the tropics, though the interpretation in 
terms of actual aerosol behavior was ambiguous as inconsis- 
tent aerosol vs. MJO-wet-phase relationships are found 
between different satellite aerosol data sets. Possible reasons 
leading to the inconsistency include the aerosol humidification 
effect, wet deposition, low-level wind effect, biological 
production, sampling effect, and cloud contamination. Tian 
et al. [2011] further investigated the MJO-related aerosol 
variability over the tropical Atlantic Ocean using the aerosol 
optical thickness (AOT) product from the MODIS (Moderate 
Resolution Imaging Spectroradiometer) on the Aqua satellite. 
They suggested that the MJO-related intraseasonal variance 
accounts for about 25% of the total AOT variance over the 
tropical Atlantic. They also found that the MJO modulates 
the Atlantic aerosol variation primarily through its influence 
on the Atlantic low-level zonal winds. Given the potential 
predictability of the MJO extending to 2-4 weeks [e.g., 
Waliser, 2005], the study by Tian et al. [2011] implies that 
the Atlantic aerosol concentration may be predictable with 
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lead times of 2—4 weeks, which in turn may lend important 
guidance to predicting air quality, dust storm activity, and 
ocean nutrient deposition over the Atlantic Ocean. 

[ 4 ] Nevertheless, Tian et al. [2011] examined only the 
total AOT anomalies and did not consider the contribution 
of different aerosol types to the total AOT anomalies. It is 
well known that the aerosol over the tropical Atlantic Ocean 
in the boreal winter season is a mixture of mineral dust from 
the Sahara desert and the Sahel region, biomass burning 
smoke from the Sahel and African savanna regions, and 
marine aerosol (primarily sea salt and secondly sulfate aerosols) 
from the ocean surface [ Kaufman et al., 2002, 2005a, 2005b], 
Because dust, biomass burning smoke and marine aerosols 
play very different roles in the radiative forcing and cloud 
formation process, it is of great interest to partition the total 
aerosol into individual aerosol components and examine the 
MJO-related variability in each aerosol type. 

[ 5 ] In this study, we use two independent satellite data sets, 
MODIS and MISR (Multi-angle Imaging SpectroRadiometer), 
to investigate the intraseasonal variability of individual aerosol 
components related to the MJO over the tropical Atlantic. 
These two data sets both have their advantages and disadvan- 
tages, thus are complementary to each other: MODIS has more 
frequent sampling than MISR due to its wider swath, but the 
scheme used to distinguish dust and smoke is subject to some 
major assumptions. MISR discriminates dust aerosol using 
actual particle shape information, however, it cannot distin- 
guish smoke from other spherical components, and it has only 
1/4 to 1/3 times sampling of MODIS. 

[6] The rest of this paper is organized as follows. Section 
2 describes the MODIS and MISR data, and the methodol- 
ogy used to derive specific aerosol components from the 
total AOT. The climatology of specific aerosol components 
as well as the comparison of dust aerosol between MODIS 
and MISR is presented in section 3 to examine the fidelity 
of the methods described in section 2. The main results of 
this paper, the MJO-related dust and smoke aerosol anoma- 
lies and correlations with MJO-related wind and precipita- 
tion anomalies, are presented in section 4. Conclusions and 
discussions are presented in section 5. 


partitioned into contributions from the dust, smoke, and 
marine aerosol components, we have the following equations: 


r Tdust ± Umoke ± X mar 

/ * T /dust X Tdi ts t T /smoke X T smoke T /mar X Ty 


( 1 ) 

( 2 ) 


[8] Rewriting (1) and (2), we get 

Tlust — [t X (/smoke /) X mar X (/smoko /mar ) | / (/smoke ./just) 

(3) 

Xsmoke | r X if /dust) X mar X (/mar /dust ) | / (/smoke /dust) 

(4) 


[ 9 ] In these equations, t denotes total midvisible AOT 
(x and AOT are used interchangeably in this paper), and the 
subscripts “dust”, “smoke”, and “mar” indicate dust, smoke, 
and marine aerosol components, respectively. With x and/ be- 
ing MODIS measurements, x dust and T smoke can be computed 
directly if/ dus t, /smoke, /man and x mar are known. 

[ 10 ] In the Kaufman et al. study, / dust , / smok e, and / mar are 
assumed to be constant, and were derived by averaging the 
MODIS/Terra/ observations over selected regions and time 
periods where one specific aerosol type dominates, with 
their uncertainties estimated from these selected measure- 
ments: / d ust = 0.5 ± 0.05, / smoke = 0.9 ± 0.05, and / mar = 0.3 
±0.1 [Kaufman et al., 2005a]. However, the actual / values 
vary with season and location [e.g., Maying et al., 2003; Yu 
et al., 2009], thus large uncertainties are expected from 
using constant / dust , / smoke , /mar when applying Kaufman’s 
formula. More discussion of this issue as well as the sensi- 
tivity of our results to the/ values adopted will be given in 
section 4.3. 

[ 11 ] The marine AOT, x mar , depends strongly on surface 
wind speed, as its primary component is sea-spray salt 
[Huang et al., 2010; Smirnov et al., 2003]; it is estimated 
using the empirical formula in Kaufman et al. [2005a] 

Xmar = 0.007 W ±0.02 (5) 


2. Methodology and Data Description 

2.1. Moderate Resolution Imaging Spectroradiometer 

[ 7 ] Previous studies [Kaufman et al., 2002, 2005a, 2005b] 
have suggested that satellite data distinguishing fine-mode 
aerosols from coarse-mode aerosols cotdd be used to sepa- 
rate the aerosol into specific types, because different aerosol 
types (e.g., smoke, dust, and marine aerosols) have different 
fine mode fraction if) values. / is the fraction of total 
midvisible AOT contributed by the fine-mode aerosols. In 
MODIS, the fine-mode aerosols refer to the aerosols with a 
size distribution of radii centered between 0. 1 and 0.25 pm, 
whereas the coarse-mode aerosols have radii centered 
between 1 and 2.5 pm. In their series of studies [Kaufman 
et al., 2002, 2005a, 2005b], they developed empirical for- 
mulae to estimate dust and smoke aerosols using MODIS/ 
Terra AOT and f observations. These formulae have since 
been widely adopted by the community to understand aerosol 
types and their climate forcing [e.g., Bellouin et al., 2005; Yu 
et al., 2009, among many others]. With the two constraints that 
both the total AOT and its fine mode fraction can be 


[ 12 ] Here W is the surface (10 m) wind speed from the 
European Centre for Medium-Range Weather Forecasts 
(ECMWF) ERA-Interim reanalysis [Dee et al., 2011], Die 
global mean value of x mar is around 0.06 ±0.005 [e.g., 
Kaufman et al., 2001]. 

[ 13 ] In this study, daily AOT and / measurements at 
0.55 pm from the MODIS/Aqua Level-3 Collection 5.1 data 
set [Remer et al., 2005; Levy et al., 2009] on 1° x 1° spatial 
grids are used. The uncertainties of x are ±(0.03 + 0.05x) 
over ocean and ±(0.05 + 0.1 5x) over land. Over ocean, the 
uncertainties of / are large for low AOT (x<0.15) but 
typically less than about 20% for large AOT (x>0.15) 
[Kleidman et al., 2005; Remer et al., 2005], Over land, 
MODIS does not provide any quantitative infonnation about the 
aerosol size [Le\y et al., 20 1 0]. As a result, we use the x and / data 
over the ocean only. The period of 4 July 2002 to 1 June 2009 is 
used for consistency with the study by Tian et al. [2011]. 

[ 14 ] Note that we use MODIS/Aqua aerosol data rather 
than MODIS/Terra data even though MISR is on board 
Terra. One reason is that we want to be consistent with the 
study in Tian et al. [201 1], More importantly, MODIS/Terra 
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and MISR aerosol retrievals have very little spatial overlap 
due to the exclusion of Sun glint regions over the ocean 
for MODIS. As a result, there are far fewer days on which 
both MODIS/Terra and MISR made observations of a given 
region, compared to the situation between MODIS/Aqua and 
MISR. Because the direct comparison between daily 
MODIS and MISR aerosol data serves as an essential part 
of this study, we chose to use MODIS/Aqua data. Neverthe- 
less, the difference between the results based on MODIS/ 
Terra and those based on MODIS/Aqua is found to be neg- 
ligible for our application. 

[is] Data rejection has been performed on the MODIS x 
and / observations as discussed in Appendix A. It is noted 
that we applied equations (3) and (4) in a stricter way than 
what was originally done by Kaufman et al. [2005a, 
2005b], Nevertheless, it should be noted that although the 
data rejection we applied removes a large number of suspi- 
cious observations, doing so does not in itself guarantee 
the accuracy of the remaining points, given that /just, /]moke, 
and f mar are assumed constant, and x mar is estimated from an 
empirical formula. For the former factor, sensitivity tests are 
performed to examine the extent to which our results are 
sensitive to the variations in f values. For the latter, the devi- 
ations of the actual x mar (mainly sea salt) compared to that 
empirically computed will be interpreted primarily as 
changes in dust amount, thus affecting the accuracy of i dust . 
However, as will be presented next, MISR can distinguish 
dust from sea salt based on retrieved particle shape, thus x dust 
derived from MISR is not subject to the ambiguity between 
dust and sea salt that affects Kaufman’s method. In this 
sense the examination of MISR aerosol observations pro- 
vides validation for the more extensive MODIS data set, in 
addition to offering actual results independent of MODIS. 

2.2. Multi-Angle Imaging SpectroRadiometer 

[ 16 ] Viewing the Earth simultaneously at nine widely 
spaced angles, MISR/Terra provides global coverage every 
7-9 days. The variations between the reflectance acquired 
from a range of observation angles can be interpreted (with 
appropriate models) in terms of aerosol properties such as 
particle size, shape, and single-scattering albedo [Kahn 
et al., 1998, 2001; Chen et al., 2008]. In particular, MISR’s 
sensitivity to the characteristics of the aerosol scattering 
phase function enables it to distinguish between the 
nonspherical and spherical particles, and thus provides a 
possible way to separate mineral dust aerosols from other 
aerosol components. A series of studies has explored the 
ability of MISR to retrieve mineral dust properties theoreti- 
cally [Kahn et al., 1997, 2001; Kalashnikova et al., 2005; 
Kalashnikova and Kahn, 2006], the sensitivity of the theo- 
retical results, as well as the application of nonspherical dust 
models for five Saharan dust field events over the Atlantic 
Ocean [Kalashnikova and Kahn, 2008], 

[ 17 ] Assuming the mineral dust is all nonspherical and 
the nonspherical part of AOT is all from dust, MISR x dust can 
be directly computed as the nonspherical fraction of the 
total AOT: 

Xdust — X X fraCnonspherical (6) 

[is] For MISR, smoke and marine aerosols (both sea salt 
and sulfate aerosols) contribute to the spherical part of the 


total AOT. (Refer to Table 3 in Kahn et al. [2001] for more 
details on shape categories of different aerosol components.) 
If T mar is again taken as computed from equation (5), x smoke 
can be approximated as the difference between the total x 
and x dust and x mar 

Xsmoke X X du st X mar (7) 

[ 19 ] However, it should be kept in mind that unlike MISR 
x dus t, the MISR x smoke derived this way, and all the MODIS 
particle type distinctions, are not independent observations 
based on actual physical constraints, and their accuracy is 
subject to the large uncertainties due to the empirically 
calculated x mar . 

[ 20 ] We use MISR-derived Level 3 daily x and 
nonspherical fraction at 0.558 pm on 0.5° x 0.5° spatial 
grids during the same time period as the MODIS data. 
MISR x and x dust are averaged onto 1° x 1° grids to com- 
pare with MODIS results. A number of validation studies 
have shown that overall, about 70% to 75% of MISR 
AOT retrievals, their retrieval errors fall within 0.05 or 
20% x AOT, and about 50% to 55% are within 0.03 or 
10% x AOT, except at sites where mixed dust and smoke 
are commonly found [Kahn et al., 2010]. Particle property 
validation suggests that expected MISR sensitivity to the 
spherical versus nonspherical particles is about 20% for 
AOT above 0.15, and diminishes for midvisible AOT 
below this value [Kahn et al., 1997; Kalashnikova et al., 
2005; Kalashnikova and Kahn, 2006], Thus, we only use 
x>0.15 to calculate x dust from equation (6). With this x 
cutoff, it is found that more than 70% of the observations 
still remain for calculation over most of the tropical 
Atlantic. Note that generally, MISR x retrievals are avail- 
able about 15% of the time, except over the convectively 
active regions (decreased to about 10%). Although it has 
never been explicitly addressed, MODIS sensitivity to 
particle properties, e.g., f also diminishes at low AOT 
(implied in Figure 7 in Kahn et al. [2009]), and MISR 
actually has greater sensitivity at low AOT than MODIS 
due to the long atmospheric paths observed by its 
steeper-viewing cameras. 

[ 21 ] Note that MISR also retrieves aerosol particle size 
information, thus conceptually it is possible to follow the 
same method utilized for MODIS to separate the total 
AOT into specific types in MISR. However MISR catego- 
rizes the aerosol particles into three bins: “fine” (particle 
radii < 0.35 pm), “medium” (radii between 0.35 and 0.7), 
and “large” (radii >0.7 pm) modes, instead of two bins as 
“fine” vs. “coarse” in MODIS. Thus, applying Kaufman’s 
fonnula to MISR would require considerable additional 
work, but not necessarily lead to greater insight, because 
similar assumptions would be required to apply the 
MODIS size-discrimination method to MISR. Furthermore, 
the different radii range for MODIS and MISR “fine” 
mode would make it impossible to cross-validate the / 
values between these two data sets. Most importantly, the 
retrieval of x dust in MISR is based on actual aerosol physical 
property, thus considered to be preferable to Kaufman’s 
method. Therefore, using MISR aerosol size information to 
separate different aerosol components following Kaufman’s 
method is not only beyond the scope of current paper, but also 
less desirable. 
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Figure 1. MODIS climatological mean (2002-2009) boreal winter (a) total x, (b) T dust , (d) x smoke , and 
(f) x mar over the tropical Atlantic Ocean, and their percentages to x ((c) x dust , (e) x smoke , and (g) x mar ). Nine- 
point spatial smoothing is applied. 


[ 22 ] Beside the MODIS and MISR aerosol data, 1° x 1° 
precipitation from Tropical Rainfall Measurement Mission 
3B42 Version 7 data set [. Huffman et al., 2007] and 1.5° 
1.5° daily horizontal winds from the ECMWF ERA- 
Interim reanalysis [Dee et al., 2011] during the same time 
period are also used. 

3. Comparison of Aerosol Components Between 
MODIS and MISR 

[ 23 ] In this section, the AOTs of individual aerosol com- 
ponents over the tropical Atlantic Ocean are examined, and 
the MODIS and MISR results are compared. Climatological 
maps are examined first, to investigate how well the methods 
presented in section 2 capture the basic features of the long- 
term mean. 


3.1. Aerosol Winter Climatology 

3.1.1. Moderate Resolution Imaging Spectroradiometer 

[ 24 ] The MODIS winter climatological mean (November 
to April during 2002-2009) aerosol maps over the tropical 
Atlantic Ocean are shown in Figure 1. The climatological 
mean x features a zonally oriented, optically thick aerosol 
plume centered at around 5°N-8°N stretching across the 
equatorial Atlantic Ocean. The magnitude and latitudinal 
extent are greatest over the eastern equatorial Atlantic with 
the maximum x up to 0.5 and gradually decrease westward 
toward the central and western equatorial Atlantic, as 
expected for an aerosol plume that originates in Africa 
(Figure la). Note that due to data rejection performed in this 
study, Figure la is highly similar but not identical to the 
color shadings in Figure lb in Tian et al. [2011] although 
they display the same quantity: the magnitude of former is 
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Figure 2. Same as Figures 2a-2e but for MISR, and only cases with MISR x>0.15 used. 


about 80% of the latter over the equatorial region. Neverthe- 
less, they have a spatial correlation of 0.98. 

[ 25 ] The spatial pattern of climatological x dust closely 
resembles the x pattern, with maximum x dust (about 0.35) 
found along the west coast of Africa, and a gradual decrease 
toward the central and western equatorial Atlantic (Figure lb). 
Over the equatorial Atlantic Ocean, where aerosol loading is 
high (x >0.15), dust is the dominant aerosol component, 
contributing more than 50%, and as much as 75%, to x 
(Figure lc). 

[ 26 ] Compared to x dust , x smoke is significantly weaker. 
Plumes of fine-mode-dominant aerosol are found originating 
from the African (biomass burning smoke) as well as the 
South American continents (smoke or air pollution) 
(Figure Id). A contribution of more than 20% to total x is 
found over the eastern tropical Atlantic, whereas the contribu- 
tion is less than 20% over the western part (Figure le). 
Figure If shows that x mar is very small (about 0.04) over the 
Atlantic intertropical convergence zone and the west coast of 
Africa because of the weak trade winds. Over the clean 
subtropical Atlantic, marine aerosol is the dominant compo- 
nent (>50%) due to a lack of dust and smoke aerosols over 
these regions. 

[ 27 ] These results indicate that the major aerosol plume over 
the equatorial Atlantic Ocean in Figure la is the dust originat- 
ing in the Sahara desert and the Sahel region, with some 
contribution from biomass burning smoke originating in the 


Sahel and African savanna. These aerosol distributions are 
generally consistent with previous observational results [e.g., 
Husar et al., 1997; Kaufman et al., 2005a; Huang et al., 
2010], lending some confidence in estimating dust and smoke 
aerosols based on equations (3) and (4). 

3.1.2. Multi-Angle Imaging SpectroRadiometer 

[ 28 ] Multi-angle Imaging SpectroRadiometer winter cli- 
matologies of x, x dust , and x smoke are shown in Figure 2. 
Recall that for these results only cases with x>0.15 are used 
due to the fact that sensitivity to shape diminishes when 
aerosol concentration is low. Different data rejection is 
applied to MODIS x, as discussed in the Appendix A, thus, 
we do not expect Figure 2a to be the same as Figure la. 
Nevertheless, we do find that they have highly consistent 
spatial patterns and comparable magnitudes. The climato- 
logical MISR x dust map (Figure 2b) is also very similar to 
the MODIS x dust (Figure lb). The MISR x and x dust are 
slightly larger than the MODIS counterparts partly due to 
MISR’s exclusion of low aerosol cases, but more important 
reasons will be addressed in section 3.2. It is seen that MISR 
Xdust contributes more than 40%, up to more than 60%, to the 
total x for cases with x>0.15, further confirming that dust is 
the dominant aerosol component over the equatorial Atlantic 
Ocean (Figure 2c). It is also noted that MISR has 10% to 
15% lower dust fraction over the equatorial Atlantic. Again, 
this difference could be partly attributed to the x cutoff 
applied to MISR, but more importantly could result from 
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a)corelation coefficient b)count 
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Figure 3. (a) Correlation coefficients between coincident MODIS x dust and MISR x dust during 2002- 
2009 over the tropical Atlantic and (b) the number of coincident days. Time series of MODIS (blue) 
and MISR (red) i dust as the function of the coincident day (c) at the point of 29.5°W, 19.5°N, averaged 
on the (d) 5° x 5° and (e) 10° x 10° boxes centered at 29. 5° W, 19.5°N. 


both uncertainties involved in the derivation of MODIS and 
MISR i dust , especially in MODIS, given the assumptions 
used. We return to this point in section 3.2. The climatolog- 
ical MISR x S moke has a 10% larger contribution to total x 
than MODIS x smo ke over the equatorial Atlantic 
(Figures 2d and 2e compared to Figures Id and le). 

[ 29 ] The overall consistency between the MODIS and 
MISR aerosol climatology, representing two independent 
satellite data sets and methods used to derive the dust and 
smoke aerosols over the tropical Atlantic Ocean, gives us 
some confidence in applying the results. However, the con- 
sistency in climatology does not guarantee their consistency 
on shorter time scales, for example, on the intraseasonal time 
scale of importance here. Therefore, in the next subsection 
we directly compare the coincident daily MODIS and MISR 
x d ust over the tropical Atlantic Ocean. 


3.2. Comparison of Coincident MODIS and MISR 
Dust Aerosols 

[ 30 ] The correlation between the daily MODIS and MISR 
Xdust for all their coincident days (minimum of 5) during 
2002-2009 is calculated over the tropical Atlantic. Here the 
term “coincident” simply means the MODIS and MISR obser- 
vations fall within a same grid box and on a same day, which is 
less strict than the usual definition used by the satellite 
community. We use the entire years for 2002-2009 instead 
of winters only to obtain as many coincident days as possible. 

[ 31 ] Overall, the MODIS and MISR x dust are well corre- 
lated in the regions of primary interest here, where dust dom- 
inates. The correlation is systematically higher in the north 
Atlantic region (typically around 0.7 or larger) than in the 
south Atlantic where the correlation is typically less than 
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a) AOT diff (MISR-MODIS) b) dust AOT diff 
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Figure 4. (a) MISR and MODIS x difference averaged over coincident days during 2002-2009. (b) Same as 
Figure 4a, but for x dust . (c) Number of coincident MODIS and MISR observations (black line), averaged 
MODIS i (thick blue line), MISR x (thick red line), MODIS x dust (thin blue line), and MISR x dust (thin red line) 
over 20°S-30°N, 60°W-20°E during 2002-2009 as a function of binned MISR x (binned by every 0.01). 


0.5, with quite a few spots less than 0.2, or even negative 
(Figure 3 a). Note that there are more coincident days over the 
north Atlantic (Figure 3b). Furthermore, the correlation is 
higher when it is closer to the African continent, and decreases 
gradually as the dust is transported away from the source re- 
gion. These results suggest that MODIS and MISR x dust agree 
with each other quite well over heavy dust load regions, 
whereas they are less consistent over the regions with less 
frequent dust occurrence or small dust aerosol concentration. 

[ 32 ] Figure 3a shows that overall, the MODIS and MISR 
x dust are highly correlated; however the correlation coefficient 
does not provide information on the x dust magnitude. Therefore, 
the time series of MODIS and MISR x dust as a function of their 
coincident days at a representative grid point (29.5°W, 
19.5°N), as well as averaged within the 5° x 5° and 10° x 10° 
grid boxes centered at this point, are shown (Figures 3c-3e). 
Again it is found that the MODIS and MISR x dust are highly 
correlated on a daily basis, and it also reveals that x dust is 
systematically greater for MISR than MODIS. 


[ 33 ] Multi-angle Imaging SpectroRadiometer x dust is 
larger than MODIS x dust over most of the tropical Atlantic 
when averaged for coincident days (about 0.03 larger aver- 
aged over the basin), except in some regions north of the 
equator (Figure 4b). This difference can be traced to the 
difference between MISR and MODIS total AOT (Figure 4a). 
Over almost the entire tropical Atlantic, coincident MISR 
AOTs are systematically larger than the MODIS ones (up to 
0.08, and about 0.04 when averaged for the Atlantic basin). 
This is consistent with previous studies in which MISR AOT 
is found to be generally larger than MODIS AOT over water 
[e.g., Abdou et al., 2005; Kahn et al., 2010]. Further examina- 
tion of coincident MODIS and MISR x and x dust binned 
against the MISR x reveals that larger MISR x is found for 
the entire x spectrum, whereas MISR x dust is larger when the 
aerosol concentration is relatively low (Figure 4c). Note that 
lower aerosol concentrations are overwhelmingly more fre- 
quent (black line in Figure 4c), thus the averaged MISR x dust 
is larger than the averaged MODIS x dust . As seen above, the 
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Figure 5. The number of strong MJO events used for MJO composite for each phase of the MJO cycle 
for (a) MOD1S and (b) MISR during the 2002-2009 boreal winters. The total number of strong MJO 
events during this period is indicated at the right upper comer of each panel. 


x d ust difference between MODIS and MISR can be traced to 
the x difference between them. However, this is not the only 
reason, as, unlike MISR x, MISR x dust is not systematically 
larger than MODIS x dust . The uncertainties involved in the der- 
ivation of both MODIS and MISR x dust inevitably contribute 
to their differences too. However it should be noted that 
although both methods have limitations, MISR dust discrimi- 
nation is based on actual retrieved particle shape information, 
whereas the accuracy of MODIS x dust is subject to the use of 
constant f values and empirically calculated x mar . Despite more 
physically robust separation of x dust from x in MISR, it has 
much less frequent sampling compared to MODIS, thus it is 
necessary to examine both data sets. Results based on these 
two complementary satellite data sets will provide more solid 
insight to the characteristics of the MJO-related dust and to 
some extent smoke variations. 

4. MJO-Related Atlantic Dust and Smoke 
AOT Anomalies 

4.1. MJO Analysis Methodology 

[ 34 ] For the MJO analysis and composite procedure, we 
use the multivariate empirical orthogonal function method 
introduced by Wheeler and Hendon [2004] and adopted 
widely by the MJO community [e.g., Waliser et al., 2009]. 
In this method, eight distinct phases of the MJO cycle are 
determined by a pair of Real-time Multivariate MJO 
(RMM) index (RMM1 and RMM2) [Wheeler and Hendon, 
2004], Typically, when enhanced MJO convection is located 
over the equatorial Indian Ocean (phase 1 and 2) and 


western hemisphere (phase 8), persistent low-level westerly 
anomalies and enhanced precipitation are found over the 
equatorial Atlantic, whereas the reverse conditions are found 
when the MJO convection is located over the maritime con- 
tinent (phase 4) and western Pacific (phase 5 and 6). 

[ 35 ] In this study, first, the intraseasonal anomalies of 
daily AOT are obtained by removing the climatological- 
mean seasonal cycle and filtering via a 30-90 day band pass 
filter. Then, a composite MJO cycle is calculated by averag- 
ing the daily anomalies that occur within each phase of the 
MJO cycle. Only days having strong MJO activity 
(RMM1~+RMM2 2 >=1) are considered for the MJO cycle 
composite. 

4.2. MJO Composites of Dust and Smoke Anomalies 

[36] Figure 5 shows the number of strong MJO days 
(RMM1 2 +RMM2 2 >=1) used for the eight-phase MJO cycle 
composite for both MODIS and MISR during 2002-2009 
boreal winters. The total number of strong MJO events dur- 
ing this period is also shown in the upper right corner of 
each panel. For both instruments, the number of strong 
MJO events used for the composite is much less than the 
actual number of total events because of satellite retrieval 
sampling issues and the data rejection applied. Generally, 
the number of MODIS events ranges from about 10 to 45, 
about as three times more than MISR, due to the much wider 
MODIS swath. 

[ 37 ] The eight-phase MJO composite maps of MODIS 
total x, x dust , and x smoke anomalies are shown in Figures 6a, 
7a, and 8a, respectively. The MISR counterparts will be 
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Figure 6. MJO composite maps of total t anomalies (multiplied by 100) for (a) MODIS and (b) MISR 
over the tropical Atlantic Ocean. Nine-point spatial smoothing is applied. 
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Figure 7. Same as Figure 6, but for x dust . 
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discussed and compared with MODIS results later. Compar- 
ing Figures 6a and 7a, MODIS t and i ( i U st have very similar 
temporal evolution and spatial patterns, except that x dust 
anomalies are slightly smaller. The spatial correlation 
between t and T dust anomalies for eight MJO phases is 
0.89. For both i and x dust , strong negative anomalies 
(as large as about —0.04) are found over the entire equatorial 
Atlantic for MJO phases 1, 2, 3, and 8. In contrast, strong 
positive anomalies (up to 0.04) are found over the equatorial 
Atlantic for MJO phases 5-6. For MJO phase 4, strong pos- 
itive anomalies occur to the north of the equator, whereas 
negative x anomalies are found to the south. The converse 
is true for the MJO phase 7. The MJO composite maps of 
total x greatly resemble those shown by Tian et ah [201 1], with 
spatial correlation 0.52 (significant at the 99.9% level), and 
have very similar magnitudes. This resemblance suggests that 
the MJO-related total AOT anomaly patterns are robust and 
not sensitive to the data sampling (fewer but more carefully 
filtered samples are used in this study). The MJO-related x smoke 
anomalies in MODIS are very weak, rarely exceeding 0.01 
(Figure 8a). The spatial correlation between x smoke and x anom- 
alies is only about 0.21. The MJO composite maps of x^, 
anomalies are not shown because x mar is linearly dependent 
on the surface wind speed, thus the MJO-related x mar anomaly 
pattern in fact reflects the wind anomalies associated with the 
MJO, which were examined in Tian et al. [201 1], Furthennore, 
it is found that the magnitudes of the x mar anomalies are 
negligible: the strongest negative/positive anomalies are about 
—0.004/0.004 (figure not shown). 

[38] The eight-phase MJO composite maps of MISR total 
x, Tdust* and T smoke anomalies are shown in Figures 6b, 7b, 


and 8b. Again, it is found that the x dust anomalies are signif- 
icantly larger than the x smoke anomalies. The x dust anomalies 
have very similar patterns to those of total x, with slightly 
smaller magnitude, whereas the x smoke anomalies are small 
and noisy. Further comparison between the MISR and 
MODIS x and x dust results (Figures 6 and 7) indicates that 
overall they exhibit very similar temporal evolution and 
anomaly patterns despite the systematically larger MISR 
anomalies compared to those of MODIS, which is likely 
the outcome of the systematically larger MISR x and x dust 
retrievals compared to MODIS over the tropical Atlantic 
Ocean, and possibly more robust dust identification for MISR, 
as discussed earlier. The rejection of all x < 0.15 data in MISR 
could also contribute to the above difference; however, MISR x 
anomalies with no data rejection do not show evident differ- 
ences (Figure not shown). Besides the magnitude difference, 
MISR x and x dust anomalies are also noisier than the MODIS 
anomalies due to lower sampling. Furthennore, although MISR 
Tsmoke anomalies are small too, they are notably stronger than 
the MODIS x smoke anomalies. Nevertheless, the consistency be- 
tween the MODIS and MISR results shown in Figures 6-8 dem- 
onstrate that dust is the dominant aerosol component on the 
intraseasonal time scale, and the MJO-related dust anomalies 
are robust, as seen from two independent sets of satellite 
observations having different strengths and limitations. 

4.3. Sensitivity of MODIS Results 

[ 39 ] The Kaufman et al. method for computing x dust and 
Tsmoke using MODIS x and f based on equations (3) and 
(4) is straightforward; however, large uncertainties in the 
computed x dust and x smoke are expected for at least several 
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Figure 9. (a) The composite MJO cycle of MODIS total x, 
x dust , Lmoke. and x mar anomalies averaged over 30°W-5°W, 
EQ-15°N for the control case, (b and c) The MJO cycle of 
MODIS idust an d T smo kc a nomalies based on the control case 
(black) and nine sensitivity test cases (color) with/n^/just, 
and of each case indicated in the legend. 


reasons: the uncertainties in the MODIS x and / retrievals, 
the uncertainties in empirically computed x mar , and most 
critically, the uncertainties resulting from assuming constant 
fnai, ,/dust? an d / moke- In this subsection, we examine the 
sensitivity of our results to the / values. We first perturb 
the f values used in this study (see section 2.1, also denoted 
as the control case) by increasing/decreasing one of them at 
a time by its uncertainty range while keeping the other two 
unchanged. These sensitivity test cases are denoted group 
1. We then test using f values derived by three other studies 
[Jones and Christopher, 2007, 2011; Yu et al., 2009], which 
are denoted group 2. These studies attempted to calibrate 
Kaufman’s technique, and they rederived the / values using 
either updated MODIS data sets, or aerosol observations 
from other satellite data sets, or using the GOCART 
(Goddard Chemistry Aerosol Radiation and Transport) 
model to locate regions where a single aerosol component 
dominates. The f values they obtained are generally consis- 
tent with what is used in the current study: f m ar , /dust? and 
/smoke are 0.25, 0.44, and 0.83, respectively from the study 
by Jones and Christopher [2007], 0.31, 0.49, and 0.78 from 
their 2011 study, and 0.45, 0.37, and 0.90 from Yu et al. 
[2009]. /values do vary spatially, and these studies find that 
they vary considerably depending on region and season. 


[ 40 ] To facilitate quantitative comparison between the 
sensitivity tests and the control case, we average the com- 
posite aerosol anomalies in each MJO phase over a represen- 
tative region (30°W-15°W, EQ-15°N) where the aerosol 
shows strong intraseasonal variations. First, the MJO cycles 
of the i, x dust , x S moke? and x mar anomalies based on the con- 
trol case (see Figures 6a, 7a, and 8a for x, x dust , and x smoke 
anomalies) are shown in Figure 9a. It is found that the mag- 
nitude of x dust anomalies averaged over this region is as large 
as that of x anomalies, whereas x smo ke and x mar have very 
small magnitudes. Sensitivity test results for x dust and x smoke 
are shown in Figures 9b and 9c, respectively. For x dLlst , over- 
all, the colored lines (nine sensitivity tests) cluster around 
the solid black curve (control case), and the MJO cycles of 
x dust anomalies based on different sensitivity tests show a co- 
herent evolution. This suggests that the MJO associated x dust 
anomalies over the tropical Atlantic are quite robust despite 
the uncertainties in using constant f values. Nevertheless, the 
lines spread in some phases. The result based on Yu et al. 
[2009] is most different from the other cases, probably be- 
cause the f values in the other studies follow the sequence 
fmar </just <Zmoke despite the deviations, whereas the order 
of/mar and/j USt is reversed in the Yu et al. [2009] study, i.e., 
the size distribution for maritime aerosol is assumed by 
Yu et al. [2009] to contain a larger AOT fraction in the 
fine mode than that for dust. The spread in group 2 is naturally 
larger than that in group 1. Furthermore, the x dust anomalies 
are more sensitive to /just? as indicated by the larger deviations 
of the dashed cyan and blue curves relative to the control 
case in phases 2, 3, and 6. 

[ 41 ] The spread of x smoke anomalies is quite large com- 
pared to their magnitude (Figure 9c). Thus, given the large 
uncertainties involved in the x smoke derivations and the small 
and noisy x smoke anomalies found in Figure 8, the numbers 
shown in Figure 9c should be more regarded as giving the 
signs of the MJO-related anomalies in each phase than 
the precise quantification of the anomalies. Nevertheless, 
the overall MJO cycle of x smoke anomalies is consistent 
among the different cases. 

4.4. Modulation of Dust and Smoke Anomalies by 
the MJO 

[ 42 ] The MJO is characterized by eastward propagation of 
enhanced or suppressed zonal wind and precipitation anom- 
alies. Over the tropical Atlantic, low-level wind and precip- 
itation are two critical parameters that lead to aerosol 
variations through their influence on aerosol emission, 
transport, and deposition. Tian et al. [2011] examined the 
lag-correlation between the total x anomalies and the low- 
troposphere zonal wind and precipitation anomalies related 
to the MJO. They found that the MJO modulates the tropical 
Atlantic aerosol primarily through the zonal wind, which can 
last as long as one MJO phase. The precipitation also 
exhibits some modulation on the aerosol, but the correlation 
is much weaker compared to the wind impact. In this sec- 
tion, we further examine the modulation of the MJO of dust 
and smoke aerosols through the lower-tropospheric zonal 
wind and precipitation. 

[ 43 ] Linear lag correlation between the MODIS x dust anom- 
alies (Figure 7a) and the low-troposphere zonal wind and pre- 
cipitation anomalies [Tian etal., 201 1, Figures 8 and 6] for the 
composite MJO cycle over the tropical Atlantic region is 
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Figure 10. Linear lag correlation coefficient between the MOD1S x dust anomalies and low-troposphere 
zonal wind anomalies (first column), x dust anomalies and precipitation anomalies (second column) for 
the composite MJO cycle over the tropical Atlantic region; similar but for i smo k e are in the right two 
columns. Negative lags means the wind or precipitation anomalies lead aerosol anomalies by a certain 
MJO phase. Only correlation coefficients above 90% confidence limit are shown. 


shown in the first two columns of Figure 10. It is found that the 
dust anomalies are strongly and negatively correlated with the 
wind anomalies over most of the tropical Atlantic region for 
lag 0, suggesting that the trade wind anomalies could modulate 
dust aerosol through westward transport, dust emission and/or 
dry deposition: trade winds produce a prevailing easterly com- 
ponent over the tropical Atlantic, thus a positive zonal wind 
anomaly means a relaxed trade wind, which conesponds to 
weakened westward dust transport, and possibly stronger dry 
deposition as well as weaker dust emission over the source re- 
gion, all resulting in negative aerosol concentration, and vice 
versa. This strong negative correlation is also evident when 
the zonal wind anomalies lead the dust anomalies by one 
MJO phase, indicating that the wind modulation lasts as long 
as one MJO phase, which is consistent with what was found 
by Tian et al. [2011]. Negative correlations between dust and 
precipitation anomalies are also found in lag 0 and lag — 1, 
but the correlations are much weaker and mostly confined 
within the convective equatorial region compared to the correla- 
tion with the wind anomalies, suggesting moderate MJO mod- 
ulation on dust aerosol through wet scavenging. Note that 
some scattered positive correlations are found especially along 
the South American coast, which might be related to the cloud 
contamination to the aerosol retrieval. 

[ 44 ] In contrast with the dust case, examination between 
the smoke anomalies and the wind and precipitation 


anomalies suggests that the MJO exhibits very little modula- 
tion of the smoke aerosol (right two columns in Figure 10). 

[ 45 ] Figure 10 shows that dust aerosol is heavily 
influenced by the MJO through both dynamical and thermo- 
dynamical processes, whereas the smoke is not. This might 
explain the relative magnitude of dust and smoke aerosol 
anomalies on the intraseasonal time scale. As shown earlier, 
the MJO composite x dust anomalies are as large as the total i 
anomalies (compare Figures 9a and 10a, and the black and 
red curves in Figure 9a). Although the MJO related i smo k e 
anomalies are noisy and small, thus difficult to measure, they 
are about 10% of the total x anomalies, as roughly inferred 
from Figure 9a. Because the above numbers are based on 
the composite of strong MJO events only, we further exam- 
ine the standard deviations of the 30-90 day filtered dust and 
smoke anomalies, together with those of the total x and x mar 
(Figure 11). It is found that the magnitude of x dust anomalies 
is comparable to or even larger than that of the total x, 
whereas the x smoke anomaly represents about 15% compared 
to that of total x, and x mar shows very little intraseasonal var- 
iability, which is overall consistent with what we found 
based on the composite of strong MJO events. Note that al- 
though the intraseasonal x dust , x smoke , and x mar anomalies add 
up to the intraseasonal total x anomaly for each day, their 
variance does not sum to the variance of total x for this case, 
suggesting the individual aerosol components are not 
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Figure 11. Standard deviation of MODIS aerosol intraseasonal anomalies (30-90 day bandpass filtered) 
during 2002-2009 boreal winter for (a) total x, (b) x dust , (d) T smoke , and (f) i mar , and their percentages to x 
((c) x dust , (e) x smoke , and (g) x mar ). Nine-point spatial smoothing is applied. 


completely independent of each other on the intraseasonal 
time scale, so the nonlinear terms do not vanish. 

[46] Recall that x dust is about 40-70% of total x, and x smoke 
and x mar contribute another around 25% each, in terms of 
winter climatological mean (Figures lc, le, and lg), which 
is quite different from the relative magnitude of x dust , x smoke , 
and x mar on the intraseasonal time scale as shown above. 
Thus, we find that although the dominance of dust over 
smoke and marine aerosols on the intraseasonal time scale 
can be inferred from their magnitudes in the climatological 
mean, the former is not necessarily proportional to the latter. 
Note that the 25% V.S. 15% difference in smoke on the 
climatological and MJO-related scales does not lead to big 
difference in absolute AOT value, actually is estimated as 
small as about 0.0015. Such a small value is likely smaller 
than the uncertainty in an individual retrieval of the smoke 


AOT. Nevertheless, it should be noted that these percentages 
are obtained by averaging all data during 2002-2009, thus 
they should be accurate enough in terms of telling the right 
direction of the change from the climatological scale to the 
MJO-related scale: For smoke, it is a decrease from 25% 
to 15%, while it is an increase from about 40-70% to about 
100% for dust. Compared to smoke, dust shows an even 
stronger fluctuation on the intraseasonal time, which is due 
to the fact that dust is more susceptible to the MJO, as it is 
more heavily influenced by both the trade wind and precip- 
itation anomalies associated with the MJO. 

5. Summary and Conclusions 

[ 47 ] Previous studies [Tian et al., 2011] found significant 
intraseasonal variability related to the MJO in the total 
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column AOT over the tropical Atlantic region. Aerosol over 
the tropical Atlantic is primarily a mixture of mineral dust, 
biomass burning smoke, and marine aerosol. Given the 
different roles in the radiative forcing and cloud formation pro- 
cesses played by these three aerosol types, as well as the poten- 
tial predictability of the MJO extending to 2—4 weeks, it is of 
great interest to further examine the MJO-related variability 
for individual aerosol types, especially dust and smoke aero- 
sols. In this study, we used two independent and complemen- 
tary satellite data sets, MODIS and MISR, to investigate the 
intraseasonal variability of individual aerosol components 
over the tropical Atlantic. 

[48] Daily MODIS/Aqua total AOT and fine-mode- 
fraction measurements are used to derive daily x dust and 
^ smoke following the method of Kaufman et al. [2005a, 
2005b]. This method contains considerable uncertainties, 
due to the assumption of constant /j ust , / moke , /mar and the 
empirical calculation of x mar . Strict data rejection has been 
applied to use Kaufman’s formula in a safer way. With 
MISR’s sensitivity to aerosol particle shape, dust and smoke 
aerosols can be distinguished using the MISR aerosol 
nonspherical fraction. Because MISR sensitivity to shape 
diminishes when aerosol concentration is low, only x >= 
0.15 data are used to compute x dust and x smoke for MISR. 
The examination of both data sets is necessary, and results 
from MODIS and MISR are complementary: MODIS pro- 
vides better spatial coverage, and MISR dust is derived from 
actual aerosol property retrieval rather than assumed aerosol- 
type-specific factors. 

[ 49 ] MODIS and MISR show a very similar dust and 
smoke winter climatologies. x dust is found to be the domi- 
nant aerosol component over the tropical Atlantic. It is larg- 
est over the eastern equatorial Atlantic (about 0.35) and 
gradually decreases toward the west. The contribution of 
T du st to total T ranges from 40% to 70%, considerably larger 
than x smoke , which contributes approximately 25%, and x mar 
contributes less than 20%. The daily MODIS and MISR x dust 
distributions are highly correlated overall, with the correla- 
tion coefficients typically about 0.7 over the North Atlantic, 
but much smaller or even negative over the South Atlantic. 
MISR x dust is found to be systematically greater than the 
coincident MODIS x dust , and this difference can be traced 
to the AOT difference between them. The consistency of 
the MODIS and MISR dust and smoke aerosol climatologies 
and daily variations lend confidence to our use of these data 
sets to investigate the relative contributions of dust and 
smoke aerosols to the total AOT variation associated with 
the MJO. However, the identification of smoke is much less 
certain than that of dust, because discrimination among fine- 
mode sea salt, sulfate, and smoke particles depends on 
assumptions for both MODIS and MISR, whereas the MISR 
dust discrimination is based on retrieved particle shape. 

[ 50 ] For MODIS, the MJO composite maps of x dust anom- 
alies are very similar to those of x anomalies, and are of 
comparable magnitude. Furthermore, the variance of the 
Xdust anomalies on the full intraseasonal time scale is found 
to be comparable or even bigger than that of the x anomalies. 
In contrast, the MJO-related x smoke anomalies are rather 
small, barely exceeding 0.01, and the x mar anomalies are 
negligible. The sensitivity study further shows that the 
MJO-related x dust and x smoke anomalies are quite robust, 
even when we perturb the / values by their uncertainty 


ranges or use different sets of f values from several indepen- 
dent studies. 

[ 51 ] Similarly, MISR also shows that the MJO composite 
maps of x dust anomalies are very similar to those of the x 
anomalies, whereas the MJO-related x smoke anomalies are 
rather small. The composite MJO cycle of x dust anomalies 
from MISR over the tropical Atlantic Ocean is consistent 
with the MODIS results, although the anomalies are much 
noisier due to less frequent sampling. The magnitude of 
MISR anomalies is again found to be systematically larger 
than that of MODIS. The MJO-related x smoke anomalies in 
MISR are overall slightly larger than for MODIS, but still 
much smaller compared to the x dust anomalies. 

[ 52 ] The consistency between the MODIS and MISR x dust 
and x smoke anomalies in terms of the evolution of the MJO 
cycle and the spatial pattern of anomalies suggests that dust 
aerosol is the dominant component on the intraseasonal time 
scale over the tropical Atlantic Ocean, and supports the other 
main conclusions of this study. 

[ 53 ] Although dust contributes about 40-70% to total x in 
terms of climatological mean, its intraseasonal variability is 
comparable to or even larger than the intraseasonal variabil- 
ity of total x. The smoke intraseasonal variability is about 
15% of that for total x, even though its climatological mean 
contributes to about 25% of the total x. This suggests that the 
MJO does not influence dust and smoke in the exactly same 
way. Examination of the lag correlation between dust and 
smoke aerosol anomalies with the low-level zonal wind 
and precipitation anomalies support this speculation: dust 
is more susceptible to the MJO as is more heavily influenced 
by the trade wind and precipitation anomaly associated with 
the MJO. The modulation of the MJO-related wind anoma- 
lies on the dust anomalies can last as long as one MJO phase. 

[ 54 ] The observational results obtained from two comple- 
mentary satellite data sets can be used to evaluate chemical 
transport models and help in model development. Further- 
more, the findings of this work have broader implications, 
related to the predictability issue. As we know, predictability 
beyond the synoptic time scale relies merely on our knowl- 
edge of the semi-periodic phenomena such as the MJO. 
Given that potential prediction of the MJO is extended to a 
few weeks currently, our finding suggests that dust has a po- 
tential prediction time scale up to a few weeks, too. This 
adds potential guidance to the prediction of phenomena 
affected by dust, such as dust storms, Atlantic tropical cyclo- 
genesis, Atlantic tropical cyclone evolution, ocean fertiliza- 
tion, and so on. 

Appendix A: MODIS Data Rejection 

[ 55 ] In this work, data rejection is performed for the MODIS 
x and / measurements for the reasons described below. 

[56] The dependence of x dust and x smoke on x and /based 
on equations (3) and (4) is illustrated in Figure Ala and 
Alb, where x mar is set to 0.06 (the approximate global mean 
value of x mar ). It is seen that the larger the x, the larger x dust 
and x smoke . Also, as / increases, x smoke increases and x dust 
decreases. It is shown that equations (3) and (4) can produce 
reasonable (nonnegative) values for both x dust and x smoke 
only when the paired x and / measurements fall within a lim- 
ited region bounded by the white dashed lines (x dust = — 0.03 
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Figure Al. Calculated (a) x dust and (b) x smoke as a function of total x and f according to equations (3) and 
(4). Here/ mar = 0.3, j /dust := 0.5, A moke = 0.9, and x mar = 0.06, as indicated by four white straight lines. The 
white dashed lines indicate x dus ,= — 0.03 and x smoke = — 0.03. Contour interval: 0.1. (c) The count of paired 
MODIS x and f measurements with respect to x and f over the tropical Atlantic (20°S-30°N, 60° W-20°E) 
for 2002-2009 boreal winters. Superimposed solid black lines are the two white dashed lines in Figures 
Ala and Alb corresponding to x dust = — 0.03 and x smoke = — 0.03. 


and x smoke = — 0.03). We relax the limit of x dust or x smoke to 
—0.03 instead of 0 because the uncertainties of the MODIS 
AOT are ±0.03 over ocean, thus a small negative AOT up 
to —0.03 is regarded as indistinguishable from the value 
0 [Reiner et al., 2005], Note that calculated x dust and x smoke 
in the range of 0 to —0.03 have been set to zero. Outside this 
region, either x dust or x smoke is too negative, and correspond- 
ingly, x smoke or x dust would be larger than the total AOT, 
which is not physical. Therefore, the x and f measurements 
giving rise to such x dust and x smoke values are rejected. In 
addition, x measurements greater than 2 are excluded, given 
the large uncertainties from possible cloud contamination 
[Zhang et al., 2005]. 

[ 57 ] The count distribution of paired MODIS x and f mea- 
surements as a function of x and f is summarized for the 
tropical Atlantic (20°S-30°N, 60°W-20°E) and for the 
2002-2009 boreal winters (October to April) in Figure Ale. 
It is found that the majority of observations fall within the 
region bounded by two lines (x dus ,= — 0.03 and x smoke = 
—0.03); however, there is also a large number of 


observations falling outside the two lines, which are rejected 
in the calculation. The rejection results from both the limita- 
tions in MODIS data and the limitations in Kaufman’s 
method. For example, in some cases x is typically between 
0. 1-0.3, which is much larger than the baseline x mar , while 
f is extremely small, nearly 0 (see bottom of Figure Ale). 
These observations are very likely artifacts of cloud contam- 
ination [e.g., Tian et al., 2008; Zhang et al., 2005; Huang 
et ah, 2012]. Possible contamination also comes from the bias 
in the sea surface reflectance [Sayer et al., 2010; Grandey 
et al., 2011], For other cases, x is also significantly larger 
than x mar and / falls between f mM and / dust , suggesting the 
aerosol is a mixture of dust and sea salt (see the area 
between the lower solid line and the f= 0.3 line in Figure Ale). 
These are likely valid observations; however, due to the 
difficulties in separating dust from sea salt over this regime 
[Kaufman et al., 2005a, 2005b], we do not include these 
observations in our calculations. 

[58] Before data rejection, MODIS aerosol observations 
are available for about 70% of the days during 2002-2009 
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boreal winters over the clear subtropical North Atlantic, and 
decreases to about 35% over the cloudy equatorial Atlantic. 
The data rejection procedure described above rejects some- 
what less than half of the observations. The most frequent 
rejections occur over the subtropical oceans, possibly due 
to the ambiguity between sea salt and dust aerosol over those 
regions. 
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